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CURRENT SYSTEM 

by 

Richard A. Goldberg 
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ABSTRACT 

From a detailed study of the momentum transport 
equation for charged gaseous fluids moving through a 
neutral gas, it is suggested that the equatorial F 
region electron density distribution (the geomagnetic 
anomaly) is simply the natural steady state distribution 
one would expect for charged fluids under the influence 
of gravitational, electric, and magnetic fields and 
production and loss, when interaction with the neutral 
medium is negligible. It is shown that this steady 
state distribution maintains itself by means of a 
longitudinal current system, and a study of the 
properties of this longitudinal current system leads 
to a technique for its complete mapping in the upper 
F region of the ionosphere. Quite importantly, it is 
shown that the latitudinal electron density distribution 
that characterizes the equatorial anomaly cannot exist 
when the longitudinal current system is absent . 

The study of the longitudinal current system also 
demonstrates the possibility of large errors in 
temperature calculations which make use of vertical 
density slope measurements to determine scale height, 
and a technique is suggested to correct such calculations. 
Finally, using cited measurements for support, it is 
possible to suggest the mechanisms responsible for the 
diurnal and asymmetric behavior of the anomaly. 
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INTRODUCTION 

The expression "geomagnetic anomaly" refers to the observed 
geomagnet ically controlled distribution of electrons in the 
equatorial F region of the ionosphere. Its properties have been 
studied by means of critical frequency, f Q F2 (Appleton, 1946; 
Rastogi, 1959; Rao , 1963; Lyon and Thomas, 1963; Rao and Malthotra, 
1964); by means of the height of the F2 layer, h ffl F2 (Thomas, 1962); 
by means of constant height profiles using ground-based sounders 
(Croom, Robbins, and Thomas, 1959 and 1960); and by means of 
Alouette Topside Sounder measurements (King et al., 1963; Lockwood 
and Nelms, 1964). The properties of this phenomenon are therefore, 
reasonably well known and will not be reviewed here. 

Many recent papers have theoretically investigated different 
aspects of the geomagnetic anomaly making use of the steady state 
continuity equation for electrons in the ionosphere with limited 
degrees of success (Kendall, 1963; Rishbeth, Lyon, Peart, 1963; 
Baxter, 1964; Bramley and Peart, 1964; Moffet and Hanson, 1965). 
Some of this work is based on ideas concerning electrodynamic 
drift presented much earlier by Martyn (1947) , and has the ultimate 
aim of completely describing the explicit causes of the geo- 
magnetic anomaly. Unfortunately, it is clear that approaches 
which attempt to investigate the anomaly in this manner are 
extremely difficult because of the many unknown quantities which 
must be included, and hence the lack of complete success at this 
time . 

On the other hand, Goldberg, Kendall and Schmerling (1964) 
and Chandra and Goldberg (1964) have studied the isothermal steady 
state transport equations for charged fluids moving in a neutral 
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medium (the advantage here being that explicit causes need not 
be known) and found that they generate, with the appropriate 
boundary conditions, accurate representations of the topside 
geomagnetic anomaly as measured by King et al., (1963). This 
semiphenomenological approach has also been successful in the 
upper bottomside F region when the isothermal restriction is 
lifted (Goldberg, 1965). 

In this work we now investigate the momentum transport 
equations in further detail to demonstrate that the geomagnetic 
anomaly appears to be nothing more than the natural steady state 
distribution of electrons under the influence of production and 
loss in the earth's magnetic and gravitational fields in a 
region where their collisions with neutral particles are relative- 
ly small in number. This interpretation also leads to plausible 
suggestions for the causes of the solstice asymmetry and the 
diurnal behavior of the anomaly. 

The analysis will demonstrate that steady state (no acceler- 
ation) requirements demand the existence of a longitudinal current 
system which not only produces the latitudinal distribution of 
the anomaly but also creates a topside vertical distribution 
different from a simple diffusive equilibrium shape. This makes 
derivations of electron temperature from the altitude profiles 
of electron density dubious in regions where latitudinal density 
gradients exist. A method will be suggested for correcting these 
measurements when information is available concerning the latitudinal 
gradient of electron density in the same region where temperature 
is being measured from vertical density slopes. 

In executing the above analysis, we will first investigate 
the fundamental derivation of the general equations used in order 
to solidify the physical interpretations and implications of the 
theory involved. This is deemed necessary since we feel that 
some of the current work in this field has lost sight of the 
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correct physical interpretations of the equations. 

THEORETICAL DEVELOPMENT 

As first proposed by Johnson (1951) and Schluter (1951), 
and later repeated by Chandra (1964) , the general equation of 
motion of a multiple component gas under the action of external 
forces may be written as 


Bv 

c 

FT 


3 v 


+ v 


Br 


s Br 


_ m, v . 

p + F + E — - sk (v -v ) 
k m s +m k k s 


( 1 ) 


where the suffixes s, k stand for particle type and 


v = macroscopic velocity 


s 


= density 


p = pressure tensor 


F = external and internal (see below) forces per unit mass 


m = mass 

' = collision frequency between the and k th particle 

The last term on the right-hand side of (1) represents collisions, 
which are assumed to be elastic and two body, and may be thought 
of as a frictional drag term. The first term on the right-hand 
side of (1) is the concentration force term, and is the exclusive 
cause of diffusive motions. 

The left-hand side of (1) is the total time derivative of 
velocity and thus represents the Newtonian acceleration of the 
gaseous fluids. Because the neglect of these terms is almost 
always assumed when studying transport problems in the ionosphere 
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with reasonable degrees of success, and because this neglect 

permits linearization of the transport equations, we will 

maintain this approach in the work which follows. 

— > 

The force term F represents all possible forces including 

S 

electrodynamic, electrostatic, gravitational, centrifugal, 

— » 

Coriolis, tidal, etc. We should note that F is not purely 

external since the electrostatic forces arise from electric fields 

of both external and internal origin. Since all conceivable 

forces are included in (1), the velocity v „ represents the entire 

th S 

macroscopic velocity of the s type particle fluid. This point 
is emphasized since it seems to be ignored in many papers dealing 
with the F region diffusion problem. 

The pressure tensor will be treated as an isotropic scalar 
by neglecting the off-diagonal viscosity elements in the tensor; 
c.f. Chandra (1964). 

Finally, in the equatorial F region of the ionosphere, we 
neglect tital, Coriolis, and centrifugal forces as a first order 
approximation, leaving us with the following equations of motion 
for electron and ion fluids passing through a neutral gas: 


n m m . n m m 

eei , \ , een„ rT, \ 

v . (v -v.) + v (v -v ) 

m +m. ei v e i n +m en e n 
e i e n 


- Vp + n m g - en (E + v X B) 
*e e e & e e 


( 2 ) 


nm.m __ n.m.m _ 

e 1 - e v . (v.-v ) + — LA-il v. (v.-v ) 
m +m. e l i e m . +m in 1 n 

e l in 


= - Vp ± + n i m i g + en. 


(E + v t X B) 


(3) 



since 
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In equations (2) and 
electrons, ions, and 
field caused by both 
electrodynamic force 
use of 


(3), the subscripts e, i, and n refer to 

— > 

neutrals, respectively; E is the electric 

— > —♦ 

external and internal sources; v X B is the 
_ s 

where B is magnetic field; and we have made 


= n 


m 


(5) 


where n is number density. It is also possible to write a 
similar equation of motion for the neutral particle flow but this 
is not needed here . 

In (2) and (3) , we now assume that 


n e - n i ' N 


and since 


m 


« m. , 

l ’ 


m 


n 


( 6 ) 


(7) 


summation of (2) and (3) gives 


N m . m 

— * — * in — * — * 

N m v (v -v ) + v. (v.-v ) 

e en v e n m.+m in i n 

l n 


= -? (p +p.)+Nm. g+JXB 
© 1 1 

where J is the current density defined by 

— > *4 —4 

J = Ne (v i -v e ) 


( 8 ) 


( 9 ) 



- 6 - 

Equation (8) assumes only one type of ionizable constituent. 

However, this is not a serious limitation since the results which 

follow can easily be extended to ionospheric regions containing 

multiple ionizable constituents. Furthermore, the F region is 

usually strongly dominated by 0 + ions and the assumption is quite 

applicable in this region at most times. 

— » — • 

We note again that the velocities v and v. in (8) are the 

e i 

total macroscopic velocities of the ion and electron fluids. 

— > 

Although the electric f ield E has been elimin ated explicitly in 
(8) , this has not remov ed its contribution to the general behavio r 

—if — ♦ 

of v^ and v.. , e.g. electrodynamic drifts. Note also that electron 
ion resistance terms have been eliminated explicitly. 

We now neglect the remaining resistance terms in (8) by 
assuming that we are sufficiently high in the ionosphere for 

and v type collisions to be negligible and thereby remove 
the effect of resistance between charged and neutral particle 
fluids. The validity of this assumption is justified at altitudes 
above the F2 maximum at the equator, and may even be justified at 
lower altitudes, c.f. Goldberg (1965). 

We can then write 


■?(p e +Pi) + Nm i g + JXBs=sO 


Next , we assume the ideal gas law 


p = n k T 
s s s 


and following Goldberg (1965) , let 


T +T . 
e i 


( 10 ) 


(ID 


( 12 ) 


H 

T 


kT 


(13) 



Then we can write (10) as 
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7N ^ t , j X B X r _ „ 

N t 2kT “ 2H“ U 

T 

"4 . _| 

where i is a unit vector in the radial direction and j is a 
velocity difference term defined by 


( 14 ) 


j 


J 

N 


e(v.-v ) 
v i e / 


(15) 


Equation (14) is the general steady state description of 

the electron and ion densities under collision-with-neutral free 

conditions and under the influence of electric, magnetic, and 

gravitational fields (by steady state, we mean dv /dt = 0). 

s 

Previously, (15) has been integrated along a field line, using 
the vertical profile of density at the equator as a boundary 
condition, to provide a general description of the topside 
equatorial geomagnetic anomaly under isothermal (Goldberg, Kendall 
and Schmerling, 1964; Chandra and Goldberg, 1964; Baxter and 
Kendall, 1965) and non- isothermal conditions (Goldberg, 1965). 

The results obtained have agreed excellently with the topside 
and f Q F2 properties of the midday equinoctial anomaly (Appleton, 
1946; Croom, Robbins, and Thomas, 1959; King et al., 1963; Lockwood 
and Nelms, 1964) and the inclusion of a realistic variable electron 
temperature has lead to an improved description of the observations 
including a partial explanation of the observed fine structure. 

The convenience of integrating along a field line is rather 
obvious. It allows us to investigate the behavior of N without 
any knowledge of j, i.e. if we take the component of (14) along 
line, we obtain 


+ 


T 



. h = 


0 


■ 'N 

N 


+ 


(16) 
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where h is a unit vector along B. Since (16) is independent of 
j , it can be integrated yielding 


N(r , 6 ) 


N ( r Q ) T ( r Q ) 

t ( r .9) 




dr 

O T 


( 17 ) 


where 0 is magnetic colatitude. Here the integration must be 
carried out along a field line that has an intersection with the 
equatorial normal at i.e., for a dipole, 

r = r csc^9 (18) 

As previously stated, evaluation of (17) requires a knowledge 
of N(r Q ) , the vertical electron density distribution at the equator, 
thus making the theory semiphenomenological . Furthermore, it 
provides us with no knowledge concerning j. In the next section, 
we shall investigate the behavior of j itself since a knowledge of 
j allows us a more useful tool in interpreting and calculating 
properties of the anomaly. 

THE CAUSES AND EFFECTS OF AN ELECTRON- ION VELOCITY DIFFERENCE, j 
In the previous section, we reviewed the properties of N 

— 4 

which can be obtained by considering a single component (h) of the 
general steady state equation (14) for an electron-ion plasma 
flowing in a neutral gas. This leaves us with two more components 
from which to obtain additional information. We now write (14) in 
the three spherical polar coordinate components as 


1 SN 
N Sr 


1ST. 1 <J x 5 >r 

t dr 2 h + 2m. gH 

T 1 T 


0 


(19) 
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1_ BN _ 1_ St_ X ^8 

Nr SS Tr 8 0 2m i gH__ 


( 20 ) 


1 8N 
Nr sin0 80 


1 S_T 

Tr sin0 80 


(j X B) # 
2m. g H 


0 


( 21 ) 


In the analysis which follows, we treat B as a pure dipole field 
having its north pole in the geographic southern hemisphere. 
Furthermore, we assume that j is sufficiently small to allow 
neglect of any influences it might have on B. Then 


B 


u M 
/ p o p 

'2TF cos 


e, 


Li, M sin 0 
° P 

In 3 ’ 

r 


0) 


( 22 ) 


where is the permanent dipole moment of B and is the magnetic 
permeability of free space. Substitution of (22) into (19), (20), 
and (21) provides 


1 BN 1 St 1 ^ 0 % sinG . 

N Sr + t Sr + 2 H t 8n H^g J 0 


(23) 


1 5N 1 St ^ 0 % cosQ 

N S8 + T se + 4n H T m i g J 0 


(24) 


+ 
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H T m i g 


u M .2, 
p sin ( 

8tt 2 


L 8 


^o M p cosesinB 

4tt 2 

r 


J 0 ] = 0 (25) 


1 M 

N S0 


1 Sjr 
T S0 
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The advantage of writing (14) in the component form given by 
(23), (24) and (25) is immediately apparent. We find that the 
special two dimensional nature of the dipole field provides us 
with two independent equations, (23) and (24), for the longitudinal 


component of ion-electron velocity difference j 


0 


This will enable 


us to obtain valuable information concerning i and hence H and 

& J 0 T 

N in the following discussions. 

For the moment, let us concentrate on (25). We restrict 
ourselves to midday (or middle of the night) periods when 31 kNt/ 50 
can be considered negligible. We then obtain 


j r = Je tan 1 < 26 ) 

where I is the dip angle of the magnetic field and we have employed 
the dip angle condition 


tan 1=2 cot 0 (27) 

We observe that unless j Q becomes infinite at the equator (a 

O 

physically unrealizable situation) , j must become zero there and 
electrons will not be able to separate from ions vertically in 
this region. This is the only information obtainable concerning 
j r and j q without further assumptions, but fortunately, neither 
of these variables appear in (23) or (24) , the expressions 
governing the distribution of N in the (r, 0) plane. 

Having exhausted the information available in (25) , let us 
return to (23) and (24) , devoting the remainder of this section 
to j . We first note that combination of these two equations 
with the elimination of j simply produces the characteristic 
form obtained with (16), thereby providing no new information. 

We also note that j is a coupling parameter between the 
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horizontal and vertical gradients of density (and temperature) . 
For the special case 


h - 0 < 28 > 

we find that the density distribution is horizontally stratified 
and has the vertical hydrostatic distribution of exponential 
decay with height, provided that we neglect gradients in temperature 
compared to gradients in density (this is quite valid in the 
topside F region) . This is possibly what one might observe above 
heights where the geomagnetic anomaly is known to occur. 

In the topside anomaly region, we observe a non-zero lati- 
tudinal distribution (^ ^-) which requires the existence of a 
j from (24) . From (23) , we see that this simultaneously alters 

'P i > VI' 

the vertical distribution (« -^) from the simple exponential decay 
form we would obtain when (28) is satisfied. We can also see this 
as follows: The solution of (23) is 


N(r, , 0 ) 


F l< 9 -^> e -J’a 2B 

i (r , 6 ,0) 


( 1 - 


^o M p sinB . v . 

4T-57I ^3- V dr 


(29) 


where is an arbitrary function and the explicit dependence of 
N on t disappears if t is a constant because the region concerned 
is isothermal. If we evaluate (29) at the equator, neglecting 0 
dependence for simplicity, we obtain 


-r r 1 a - ° P J «)dr 

N(r,n/2) = e Ja ' 2H 7 4n "i® ^3 (30) 

Thus, the vertical distribution at the equator is the simple 
exponential decay law perturbed by a term containing j . In the 



12 


region of the F2 peak, j must have sufficient influence to distort 
the distribution into the required shape. 

The physics involved can be thought of as follows: Since 

the steady state transport equations describe the geomagnetic 
anomaly, the anomaly is simply the natural steady state distribution 
of the electron density to be expected in a gravitational field 
under the influence of a dipole magnetic field and the effects 
of production and loss. Production and loss influence the 
distribution by contributing to the existence of the F2 peak. 

This in conjunction with the geomagnetic field requires a 
longitudinal current system, necessary to preserve the vertical 
distribution forced by production, loss, and transport combined. 
However, the existence of a longitudinal current system in the 
geomagnetic field simultaneously requires the latitudinal 
distribution of electron density known as the anomaly. The 
latitudinal peaks of the anomaly have special significance also, 
as will become evident from the analysis in the next section. 

Experimental evidence indicating the stability of the anomaly 
has been given in studies of f Q F2 by (Rastogi , 1959; Lyon and 
Thomas, 1963; Rao and Malthotra, 1964). We find that the anomaly 
forms at about 8:00 to 10:00 am local time, depending on meridian, 
and then preserves its shape and amplitude throughout the day. 

After dark, the peaks remain relatively stable in latitudinal 
position but decay in amplitude. As morning progresses, the 
decay reaches completion and the sun begins to appear. This 
creates a non-steady state distribution which masks any effects 
associated with the anomaly. Finally, as transients disappear 
and the steady state reappears, the anomaly once again forms as 
the natural steady state distribution. Such empirical studies 
demonstrate that the anomaly is, as indicated here theoretically, 
the natural distribution of electrons and ions, one which 
disappears only when transients in the distribution caused by 
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sunrise temporarily remove steady state conditions and mask the 
visible anomaly. 

We have shown at this point that the density in the (r,0) 

plane cannot be affected by currents j and j. but only by j 

r y 0 

(see (23), (24), and (25)). Hence, studies which neglect 

the perpendicular component of velocity in the continuity equation 

are incapable of generating the anomaly. 

From (23) , (24) , and (29) we note that the vertical distri- 
bution of electron density is not influenced by scale height (and 
thus temperature) alone. We find that the simple scale height 

slope is distorted by a term containing i.. A simple order of 

0 

magnitude calculation shows that a velocity separation of electrons 
and ions (v.-v ) of only 3 or 4 cm/sec is sufficient to either 

cancel (j positive) or be equal to (j , negative) the effect of 

0 0 

gravity. In a region where j is positive, i.e. between the 
latitudinal peaks of the anomaly and the equator, one will measure 
a vertical slope leading to a scale height larger than the true 
scale height. Thus, the measured temperatures will be larger 
than the true temperatures. At midlatitudes (higher than the 
position of the latitudinal peaks) or at heights above the anomaly 
where we observe a steady decay of density with latitude, is 
negative and one will measure effective scale heights leading 
to deduced temperatures smaller than the true temperatures. 


The magnitude of j can be obtained by measuring horizontal 

1 ^ XT r 

slopes of 4 and neglecting 9 dependence of t. Using (24) 
such calculations provide results in the cm/sec order of magnitude 
range further demonstrating that is not an insignificant 
quantity . In the next section, we will show how it is possible 
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thereby allowing us to correct vertical slope measurements for 
the effects of j and thus obtain true scale heights and 
temperatures . 

THE PROPERTIES OF i , 

0 

At this stage of development, no mention has been made of 

the properties of i which must exist if (23) and (24) are to 

0 

be satisfied consistently. We shall find that i is not 

0 

completely arbitrary, but instead possesses properties extremely 
useful in separating the true scale height from the measured 
vertical slope, thereby allowing us to obtain more accurate 
determinations of the electron temperature in the topside 
ionosphere . 

The consistency requirement for equations (23) and (24) is 
easily obtained by rewriting them as 


dlnNT = ^o M p sine ^ 1 

9r 8tt kT J 0 " '2IT 


SlnNT 

Te 


^o M p cos9 • 

Irr fer 2~ ^0 
r ^ 


Then, 


d r% M p sin0 
S0 L8 tt kT 3~ J 0 


-^- 1 - 


' (InNT ) _ S f ^o^p cos 8 . 

Sr 50 9 r *- 4tt kT 2 J 0' 


In the equatorial topside, the medium is to a first approxi- 
mation isothermal, and (33) becomes 


d j 


0 


d J 


d 0 


+ 2r = 3 j 

Sr J 0 


(31) 


(32) 


(33) 


tan 8 


(34) 
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where we have also restricted ourselves to a region of constant 
hk so that 

3H 

59 1 « 0 (35) 

Equation (34) represents the constraint equation on j . It 
can be solved by separation of variables yielding 

j 0 = f K G (X) r (3 ~ x)/2 sin X 6 d \ (36) 

where \ is the separation constant and G(X) is an arbitrary 
function. In order to guarantee j finite as r - ® , (36) 
requires X^3. Then, with the following substitutions, 



TD 

II 

1 

CO 

(37) 


Y (0) = G (8+3) 

(38) 

and 

i 



u = ln(r 2 /sin0) 

(39) 

we have 


j = sin e i n Y (8) e dp 

0 O 

(40) 


Thus, since the integral in (40) is a Laplace transform, j 
is a rather arbitrary function of u depending on Y(0)„ 

Let us now consider lines of constant u, which are clearly 
the dipole field lines specified by (18). Since u is constant 
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along these paths, we can immediately conclude that 

3 0 = C sin 3 0 = C(r/r o ) 3/2 (41) . 

along a field line, where C is a constant. By specifying j 
at some given point on a field line, we can then obtain j and 
d . -k^NT _ at all other points on the same field line using (41) 
and (32) . 

The above discussion then permits the following general 
procedure. If N and t are known for a single constant height 
profile, we can calculate j along this profile from the 

logarithmic slope of Nt . Then j can be evaluated in the 

Q SlnNT 

entire (r,0) plane using (41). Finally, at any point that — 

is known from vertical measurements, we can evaluate the true 

scale height and temperature using (31). Naturally, great care 

must be taken in the evaluation of results using such a technique 

due to the large errors which can result from inaccurate slope 

determinations. Furthermore, because (41) depends on the 

assumption of isothermallity , this technique is restricted to 

regions where logarithmic gradients of N are much larger than 

those of t . 

Finally, (32) and (41) indicate the physical effect occurring 
at the field line along which the angular peaks lie, viz. that 
along this field line, j is zero. Below this latitude is 
positive and hence west to east; above this latitude, j is 
negative and hence east to west in direction. 

SOME FURTHER PROPERTIES OF THE GEOMAGNETIC ANOMALY 

Because the anomaly can be thought of as the natural steady 
state distribution of electrons in the vicinity of the geomagnetic 
equator, it is possible to explain another of its properties. 

Farley (1965), using incoherent backscatter data, has observed 
that the nighttime peak of the F2 layer often maintains itself 


17 


at the equator until local midnight at which time it begins to drop 
into the higher loss regions. Studies of f Q F2 by Rastogi (1959), 
Lyon and Thomas Q.963), Rao and Malthotra (L964), and the author 
(in process) indicate that the latitudinal distribution maintains 
itself till approximately the same time period, at which time the 
angular peaks begin to decay. 

The time constant of loss appears relatively long due to 
the length of time before observable decay begins to take place. 

It therefore appears quite plausible to think of the density 
distribution as a very slowly varying function in time always 
extremely close to steady state. Under this assumption, the 
mechanism causing the layer drop at the equator becomes evident . 
The measurements of Thomas (1962) and the theory of Goldberg, 
Kendall and Schmerling (1964) and Goldberg (1965) indicate that 
the steady state vertical F2 peak be much higher at the equator 
than at latitudes where the latitudinal f Q F2 peaks are observed. 

Since the latitudinal peaks of f Q F2 are in a higher loss 
rate region than the equatorial minimum they decay faster and 
cause the steady state distribution to slowly readjust to 
compensate for the fadeout of these latitudinal peaks. This 
requires a lowering of the F2 layer at the equator, c.f. Goldberg, 
Kendall, and Schmerling (1964) and Goldberg (1965). Thus, the 
drop of the layer at the equator is nothing more than a slowly 
varying readjustment of the entire F region distribution, never 
deviating much from the steady state, to compensate for the 
more rapid initial decay of the F2 layer at latitudes where the 
latitudinal peaks in the anomaly occur than at the equator. 
Naturally, as the equatorial layer finally lowers to more rapid 
loss rate heights equivalent to the height of the layer at higher 
latitudes, we would expect a more rapid decay at the equator, 
and this is what Farley (1965) has observed. 

The geomagnetic anomaly is also known to possess asymmetric 
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properties which are most exaggerated during months of solstice 
(Croom, Robbins, and Thomas, 1960; Lyon and Thomas, 1963), 

During such a period, the asymmetry exhibits itself as a higher 
density peak in one hemisphere. This behavior is accounted for 
in (17) by means of T(r,0), which will certainly exhibit asymmetric 
properties when heating of one hemisphere is larger than the other. 
Although (17) does not explain or define the behavior of T(r,0), 
it allows us to calculate it as a function of 0. It can be shown 
that N(r,0) is far more sensitive to fluctuations of T(r,0) 
exhibited in the coefficient of the exponential than in H . 

This allows us to think of N(r,0) as being inversely proportional 
to T(r,0) to first order. It is then a simple matter to take a 
constant height profile of the asymmetric anomaly and calculate 
the temperature distribution necessary to make the distribution 
deviate from the symmetric case. Naturally, this requires an 
approximation of the distribution which would be present during 
equinoctial symmetry, but it does not seem unreasonable to 
select the average of the asymmetric values in the two hemispheres 
as a first approximation. 

SUMMARY AND CONCLUSIONS 

We have suggested that the geomagnetic anomaly is nothing 
more than the natural steady state distribution of electrons under 
the influence of production-loss effects in the earth's magnetic 
and gravitational fields and in a nearly collision free (charged- 
neutral particle type) medium. Furthermore, properties of both 
the vertical and latitudinal distribution satisfy this explanation. 
Experimental evidence has been cited to justify these conclusions. 

From this analysis, we have also been able to offer plausible 
explanations for the asymmetric properties observed during solstice 
and the nighttime behavior of the F2 layer observed at the equator. 

We have also suggested a method for evaluating latitudinal temperature 
variations . 
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The above analysis has been obtained by studying the 
momentum transport equations as applied to the equatorial 
ionosphere in greater detail than any time previously. This 
has lead to information concerning the current density system 
which must exist to maintain both the observed vertical and 
latitudinal electron density distribution. The analysis has 
displayed the character of magnetic influence on the density 
distribution, by proving that the (r,0) plane density 
distribution is quite independent of currents in the (r,0) plane 
but depends only on longitudinal currents . 

The study of j has also shown that any region possessing 
a latitudinal gradient in density must simultaneously have its 
vertical distribution distorted from the pure diffusive equilibrium 
type shape. This could lead to rather large errors in temperatures 
obtained using vertical slope techniques in the topside F region. 
However, it is possible to estimate the errors involved, and a 
method to evaluate such errors has been suggested. 

The fact that the electron density distribution in the (r,0) 
plane is influenced directly by variations in the longitudinal 
component of current lends itself nicely to the possibility of 
explaining more local type phenomena (such as Spread F) . This 
could be achieved by study of various mechanisms capable of 
causing perturbations in the longitudinal current system thereby 
altering N(r , 0) . 
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